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Abstract. The Urban Areas Working Group is part of the International Atomic 
Energy Agency's EMRAS II (Environmental Modelling for Radiation Safety) 
Programme. The goal of this Working Group is to test and improve the capabilities 
of models used in assessment of radioactive contamination in urban settings, 
including dispersion and deposition events, short- and long-term contaminant 


redistribution following deposition events, and potential countermeasures or 
remediation efforts for reducing human exposures and doses. The Working Group 
has developed three modeling exercises, which are designed to permit 
intercomparison of model predictions and, in one case, comparison of model 
predictions with measurements. This paper describes the scenarios and provides 
comparisons of initial modeling results. Reasons for similarities and discrepancies 
among model predictions are discussed in terms of the modeling approaches, 
models, and parameter values used by different assessors. Preliminary conclusions 
emphasize the value of explaining individual approaches and the importance of 
understanding the effects of different assumptions and parameter values on the 
modeling results. 


1. INTRODUCTION 


The Urban Areas Working Group was organized within the International Atomic 
Energy Agency's EMRAS II (Environmental Modelling for Radiation Safety) 
Programme, as part of a theme entitled "Approaches for Assessing Emergency 
Situations." The Working Group is building on the work done by the Urban 
Remediation Working Group of the first phase of the EMRAS Programme [1-5]. 
The goal of the Urban Areas Working Group is to test and improve the 
capabilities of models used in assessment of radioactive contamination in urban 
settings, including dispersion and deposition events, short- and long-term 
contaminant redistribution following deposition events, and potential 
countermeasures or remediation efforts for reducing human exposures and 
doses. 

The Urban Areas Working Group has developed three modeling exercises, 
which are designed to permit intercomparison of model predictions. In one 
case, comparison of model predictions with measurements is also possible for 
several endpoints. Reasons for similarities and discrepancies among model 
predictions are discussed in terms of the modeling approaches, models, and 
parameter values used by different assessors. An important objective is the 
identification of areas in which models or selection of parameter values could be 
improved. 


2. SHORT-RANGE ATMOSPHERIC DISPERSION EXERCISE 


The short-range atmospheric dispersion exercise is based on data from several 
field tests performed by the National Radiation Protection Institute on a test 
area belonging to the National Institute for Nuclear, Chemical and Biological 
Protection in Kamenna, near Prague, Czech Republic [6]. The exercise permits 
comparison of model predictions with measurements of surface contamination, 
time-integrated air concentrations, and dose rates, up to 50 m downwind. 
Intercomparisons of model predictions are used for distances up to 2000 m 
downwind and for additional modeling endpoints. 

In these field tests, a short-lived radionuclide (Tc-99m) in liquid form was 
spread by detonation of a small amount of explosive in an open field (flat 
terrain) or in an open field with some simulated structures. Measurements 
included dose rates, surface contamination, air concentrations, particle size 
distributions, time-distributions of dust particles in air, and thermo-camera 
snapshots. The test area was selected for a stable wind direction under usual 
meteorological conditions. 

Four individual field tests were used in the exercise; selected initial 
conditions for these tests are summarized in Table 1. Participants in the 
modeling exercise were asked to submit predictions for surface contamination 
and dose rates as a function of distance, air concentrations as a function of 
height and distance. Participants were provided with all available 
measurements for the first two tests, to permit an opportunity for calibration of 
models if desired. For the third and fourth tests, participants were asked to 
submit model predictions before having access to measurements of the 
modeling endpoints. 


Table 1. Summary information for the field tests. 


Test 1 Test 2 Test 3 Test 4 


Date 6 December 15 May 2008 5 May 2009 14 July 2009 


2007 


(blowing from) 


Activity (MBq) 780 1058 1222 1088 
Volume of LSL 6 mL 6 mL 6 mL 
activity 

Type and | Permon 10T, 30 | Vesuvit TN, 20 Permon 10T, Permon 10T, 
amount of g g 350g 350 g 
explosive used 

Temperature 6.7-7.1 22.2-22.3 9.7-10.4 25.4-25.6 
(°C) 

Wind speed at 7.9-16.2 1.2-6.6 1.4-16.2 0-3.2 

2 m height 

(km/h) 

Wind direction SW-SSW S-SSW SW-NNE WSW-SSE 


Eight participants have submitted preliminary calculations for the short- 
range exercise. Table 2 summarizes selected information about the models and 
their parameterization. The models represented three main types of approaches 
to modeling atmospheric dispersion and had been developed for a variety of 
purposes. Although there was much discussion about standardizing some of the 
input information, participants generally used their own interpretations of the 
information. Examples of predicted deposition down the center line of the grid 


or plume are shown in Fig. 1 for Tests 3 and 4. 


Results for Test 3 were 


generally closer to each other and to the measurements than those for Test 4, 
which had less stable wind conditions and downwind obstacles. 


Table 2. Comparison of models and selected parameters used in the short-range 
atmospheric dispersion exercise. 
Model Type of model | Stability classes Wind speed Dry deposition 
(Participant) (m/s) velocity (m/s) 
ADDAM/CSA- Gaussian Test 3: Class C Test 3: 2.7 1E-1 
ERM (Chouhan) Test 4: Class A Test 4: 0.726 
Hotspot 2.07.1 Gaussian Test 3: Class D Test 3: 1.5 Respirable 
(Charnock) Test 4: Class C Test 4: 0.4 fraction, 1E-4; 
nonrespirable 
fraction, 4E-1 
Hotspot Gaussian Test 2: Class B Test 2: 0.6 Respirable 
(Trifunovic) Test 3: Class D Test 3: 1.3 fraction, 8E-4 
Test 4: Class C Test: 0.1 
RDD MMC Lagrangian Test 1: Class C Test 1: 4.00 0.2 pm, 5.0E-3 
(Duran) Test 2: Class A | Test 2: 0.59 1.0 pm, 1.5E-4 
Test 3: Class B Test 3: 1.30 8.0 pm, 1.0E-3 
Test 4: Class A Test 4: 0.20 20.0 pm, 8.0E-3 
University of Lagrangian Not applicable Time- Not applicable 
Seville dependent 
(Periáñez) measurements 
LASAIR Lagrangian Test 1: Class D Test 1: 0-6.3 < 0.39 pm, 5E- 
(Walter) Test 2: Class Test 2: 0.28- 5 
B-C 1.85 0.39-1.3 um, 
Test 3: Class D | Test 3: 0.9-7.2 1.5E-5 
Test 4: Class C Test: 0-4.9 1.3-10.2 um, 
1E-3 
> 10.2 pm, 8E- 
3 
CFD (de With) Computational Not applicable Steady-state 0.2 um, 5.0E-5 
fluid dynamics conditions 1.0 um, 1.5E-4 


(mean values) 8.0 um, 1.0E-3 
20.0 um, 8.0E-3 
CLMM (Fuka) Atmospheric Not applicable Test 1: 3.8 Not applicable 
computational Test 2: 0.77 
fluid dynamics Test 3: 2.3 
Test: 0.4 
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Figure 1. Examples of model predictions and measured deposition down the center line 
of the grid (x = 0; ADDAM, CFD, CLMM, USEV, measurements) or the plume center line 
(Hotspot) for Tests 3 and 4. 


3. MID-RANGE ATMOSPHFRIC DISPERSION EXERCISE 


The mid-range atmospheric dispersion exercise is based on a hypothetical 
accident at a nuclear power plant and the resulting predicted deposition in 
urban environments up to 70 km downwind. The scenario assumes a 1-hour 
release from a rupture of a steam generator tube (based on an accident scenario 
developed by ISRN) and uses actual geographic and meteorological information 
for the Trillo nuclear power plant in central Spain. The exercise considers two 
sets of atmospheric stability conditions, stable (Class E) and neutral (Class D), 
and uses wind fields 10 m above the ground. The radionuclides are assumed to 
be released as gas, with a release height of 50 m. Only dry deposition is 
considered. Although a variety of radionuclides would be released in such a 
situation, for modeling purposes only two (Cs-137 and 1-131) are considered. 
Time-dependent release rates were provided. 

Participants were asked to carry out a simulation for a 10-h period and to 
provide estimates of deposited activity at the end of the period, time-integrated 
activity concentrations in air, and time series of activity concentrations in air at 
selected locations. This is a model intercomparison exercise for all endpoints. 

Five participants have submitted preliminary calculations for the mid-range 
exercise. Table 3 summarizes selected information about the models and their 
parameterization. The models represented two main types of approaches to 
modeling atmospheric dispersion. Examples of predicted deposition of Cs-137 
at two downwind locations are shown in Fig. 2. JRODOS was used only for the 
stable atmospheric conditions. RASCAL's predictions did not extend as far as 
Madrid. Under neutral atmospheric conditions, plumes predicted by ADDAM 
and the University of Seville model (not shown) missed or barely intersected 
Madrid. 


Table 3. Comparison of models and parameters used in the mid-range atmospheric 
dispersion exercise. 


Model Type of model Release Wind speed Dry deposition 
(Participant) time step velocity (m/s) 
ADDAM Gaussian 1h wind vectors Cs-137, 0.01 
(Chouhan) summed outside 1-131, 0.008 
the code 
Hotspot Gaussian 1h Class E: 3.0 m/s Cs-137, 0.0004 
(Trifunović) Class D: 6.0 m/s I-131, 0.0022 
JRODOS Gaussian + 30 min wind fields as calculated 
(Sdouz) simplified puff provided internally 
RASCAL 3.0.3 Gaussian plume 1h limited number of Cs-137, 0.003 
(Mancini) and Lagrangian wind vectors used 1-131, 0.003 
puff 
University of Lagrangian 1 min wind fields as not applicable 
Seville provided 
(Periáñez) 
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Figure 2. Examples of model predictions for deposition of Cs-137 at selected locations 
downwind from the Trillo NPP. Guadalajara is 46 km from Trillo, and Madrid is 70 km 
away. 


4. CONTAMINANT TRANSPORT AND COUNTERMEASURES EXERCISE 


The contaminant transport and countermeasures exercise starts with an 
assumed concentration of Co-60 or Pu-239 in air, in parts of a city (Seoul) for 
which detailed geographic and building information is available. Participants 
were asked to predict deposition for several kinds of weather conditions (dry, 
light rain, and heavy rain) and for both a business area (buildings and asphalt) 
and a park area. Additional endpoints for model intercomparison include 
contamination densities, dose rates, countermeasure effectiveness, and doses 
for specified reference individuals. 

Three participants have submitted preliminary calculations for the 
countermeasures exercise. Table 4 summarizes selected information about the 
models and their use in this exercise. Both METRO-K and CPHR started with 
the air concentration as provided in the scenario description. ERMIN started 
with the deposition on a lawn as predicted by METRO-K for a given set of 
conditions. Fig. 3 shows the predicted dose rate from Co-60 for an outdoor 
location in the business area, for a contamination event occurring under dry 
conditions in the summer, together with the % contributions to the dose rate 
from the most important surfaces. Fig. 4 shows the predicted external doses 
over 1 year or 5 years from a defined situation including that same outdoor 
location and nearby indoor locations. Fig. 4 also shows the % reduction in dose 
expected for several specified countermeasures. 

Some differences in the results are attributable to differences in how the 
models handled different parts of the situation. For example, roads, trees, and 
soil/grass had different importance in the different models (Fig. 3); therefore, 
the models differed in the dose reduction achieved by countermeasures 
involving these surfaces (Fig. 4). The one model (CPHR) that included the 
contribution to external dose of the plume also predicted the highest dose 
reduction due to relocation. 


Table 4. Comparison of models and parameters used in the contaminant transport and 
countermeasures exercise. 


Model Deposition Weathering Indoor Trees 
(Participa contamination 
nt) 
METRO- | dry deposition by | short-term and not included deciduous; date 
K surface; wet long-term for leaf fall not 
(Hwang) | deposition from removal rates by specified; 
daily rainfall and | surface deposition on 
washout ratio, trees in winter 
with retained 1096 of 
fraction deposition in 
summer 
ERMIN deposition on surface-specific included, from deciduous and 
(Charnoc | lawn from empirical penetration of coniferous; 
k) METRO-K; other | weathering building; simple specified date for 
surfaces relative | functions; empirical leaf fall 
to lawn, wet or movement down | retention 
dry soil column function for 
generic indoor 
surfaces 
CPHR dry deposition half-lives indoor air but not | deciduous, no 
(Tomas) velocity; washout | depending on surfaces; contamination 
coefficient for surface filtration factor after first leaf fall 


| wet deposition 


5. SUMMARY AND CONCLUSIONS 


For all test exercises, more detailed analysis is in progress to explain the 
similarities and differences among models and between model predictions and 
measurements. For each test exercise, participants have started with the same 
information, and in several cases there was considerable discussion and 
attempts to agree on specific input information, such as the dimensions of the 
initial plume for the field tests. However, differences in individual approaches 
and models remain. Areas of differences include the type of dispersion 
modeling incorporated into a code, the type of input information required by a 
model, components included in the model (e.g. indoor surfaces in the 
countermeasures models), and values selected for various parameters (e.g., size 
distribution of particles, deposition velocities, weathering coefficients). To 
understand either similarities or differences among modeling results, it is 
essential to explain individual approaches and to understand the effects of 
different assumptions and parameter values on the modeling results. 


Hl paved E wall 
Location 4, Co-60, June, d| E Ex... O a 


m 
Q 


Cobalt-60, dry, summer, Location 


E o Year 0 METRO-K 
~ l —9— METRO-K| F 
è & --Er- ERMIN Year 1 
E agit. ©- CPHR Year 5 
d R E 4 
Ss J 
: Year 0 es |" 
ELE 4 Yers] a — 
i E 4 
$ 
* Year 0 dran 
3 Year 1|. | 
10 i Vs i E 
o - N o + roy o Year 5 
a E. E. E. E E. E. 
2 2 2 2 2 2 2 0 20 40 60 80 100 


% contribution to dose rate 
Figure 3. Predicted external dose rate (left) and % contribution to dose rate of various 
surfaces (right) from Co-60 at an outdoor location in a business area, for a deposition 
event occurring under dry conditions during summer. 
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Figure 4. Predicted external dose (for a defined exposure situation) from Co-60 over 1 
year or 5 years (left) assuming no action (no countermeasures), and % reduction in 
predicted external dose expected over 1 year (right) assuming application of individual 
countermeasures. 
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